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Limits on the effective resolution of many optical near-field experiments 
are investigated. The results are applicable to variants of total internal reflec- 
tion microscopy (TIRM), photon scanning optical microscopy (PSTM) and 
near-field optical scanning microscopy (NSOM) in which the sample is weakly 
scattering and the direction of illumination may be controlled. Analytical ex- 
pressions for the variance of the estimate of the complex susceptibility of an 
unknown two-dimensional object as a function of spatial frequency are ob- 
tained for Gaussian and Poisson noise models and a model-independent mea- 
sure is examined. The results are used to explore the transition from near-zone 
to far-zone detection. It is demonstrated that the information content of the 
measurements made at a distance of even one wavelength away from the sam- 
ple is already not much different from the information content of the far-field. 

(c) 2004 Optical Society of America 

OCIS codes: 110.0180, 180.5810, 290.3200, 070.6020, 110.3000. 

1. Introduction 

Near-field microscopy 1,2 has emerged as an increasingly important imaging technol- 
ogy in the last two decades. The principal advantage of near-ficld optical methods 
is the ability to provide images with resolution on scales much smaller than the 
Rayleigh- Abbe resolution limit. 3 Synge 4 first suggested the basic method in which 
a thin sample is illuminated through a subwavelength aperture that is scanned very 
close (less than a wavelength away) to the sample while the amount of transmitted 
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light is recorded as a function of position. This method is now known as near- field 
scanning optical microscopy (NSOM ), 5-8 and is often practiced by illuminating the 
sample and collecting the light: with a subwavelength aperture at, the end of a tapered 
optical fiber. Photon scanning tunneling microscopy (PSTM) is a related method in 
which the sample is illuminated by an evanescent wave created by total internal re- 
flection at a prism face and the subsequently scattered light is collected by a tapered 
fiber as in NSOM. 

Total internal reflection microscopy (TIRM) also takes advantage of sample illu- 
mination by evanescent waves at a prism face. Unlike NSOM or PSTM, in TIRM the 
scattered radiation is measured in the far-field. TIRM has been used primarily as a 
surface inspection technique , 9,10 however, recently there has been interest in taking 
advantage of the high spatial frequency content of the evanescent wave 1 1 to perform 
direct imaging with transverse resolution surpassing the Rayleigh- Abbe limit . 12 

Quantitative interpretation of NSOM, PSTM and TIRM images is sometimes 
difficult. Under certain simplifying assumptions such as homogeneity of the mate- 
rial dielectric susceptibility . 6 the measurements may be related to topography. Under 
more general conditions, a solution of the linearized scattering problem 13-16 may be 
employed to relate near-field optical measurements to object structure and composi- 
tion. An important question remains: What is the resolution that a given experiment 
can hope to achieve? 

This paper addresses the question of resolution through an analysis of the vari- 
ance in the Fourier components of the estimate of the sample. The linearized inverse 
scattering solution 13-16 is applied to the case of a thin, weakly scattering sample and 
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an analytic expression for the error covariance matrix of the estimate of the sam- 
ple susceptibility is obtained for Gaussian and Poisson noise models. Furthermore, 
a noise-model independent bound on the estimate error is obtained. These formulae 
are then applied to PSTM for increasing distances between the sample and measure- 
ment plane until the experiment is essentially TIRM. It is seen that in the ncar-zonc 
one may surpass the usual diffraction limit with the bounds being set by the signal- 
to-noisc ratio. It is shown that the information content in the plane that is only a 
wavelength away from the sample is already quite similar to that of the far-held, i.e., 
one obtains the usual Ewald sphere of reflection. 3 

2. Scattering model 

Consider an experiment in which a monochromatic field is incident on a dielectric 
medium with susceptibility //(r). The field incident on the sample is taken to be a 
plane wave. One half space, taken to be z > 0, has the vacuum index of refraction 
while the z < 0 half space has an index of refraction n. Only nonmagnetic materials 
are considered and so it is sufficient to treat only the electric held E. The field satisfies 
the equation 

V x V x E(r) — n 2 (z)k% E(r) = 47r/c§7?(r)E(r), (1) 

it 

where ka is the free space wave number and n(z) is the z dependent background index 
of refraction as decribed above. 

The sample is assumed to be of constant thickness A z % where A z is much less 
than a wavelength, and to depend only on the transverse spatial coordinate p, where 
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r = (p, z) so 


77 (r) = 7]{p) for 0 < z < Az 
0 for 0 > Az . 

The field is taken to consist of two parts, E = E' + E^ , where E l is the incident 
field satisfying Eq. (1) with 77 ( r) = 0. The scattered field E 3 arises as a result of 
77 (r) / 0. The incident field is a (possibly evanescent) plane wave 

E*(r) = e*exp [xk(q*) ■ r] , (2) 

where q* is the trail verse wave vector of the incident wave and e* is the complex field 
vector such that |e *( 2 has units of energy/ (time • area). The wave vector k is specified 
by the transverse components q such that 

k(q) = (q,* 2 (q)) , and (3) 

q) = \Ao - IqI 2 • (4) 

The modes for which |q| < k 0 are liornogeneous, or propagating, plane waves. When 
|q| > k 0 the plane wave is evanescent, decaying exponentially with increasing values 
of 2 . These waves are super oscillatory in the transverse plane and thus provide a, 
means to probe the high spatial frequency structure of the sample. Keeping only the 
first term in the Born series , 3 the scattered field is given by the expression : 1 

K( r ) = k o [ dV Gaf}( r, r )e' 0 exp [iq< • p ' T k z (q i)z'\ r;(r') ? (5) 

where G is the 1 /2-space Green's tensor , 1 and p' is equal to the transverse components 
of r'. 
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The measurements are made in the plane 2 = z d , which is illuminated by a 
reference wave with polarization e r 

E r (r) = e r exp [ik(q r ) • r] . (6) 

Note that in the NSOM and PSTM imaging modalities, the incident field also plays 
the role of the reference field. The total field in the detector plane is given by the 
superposition of the reference and scattered fields, and the resulting intensity, called 
the iolal intensity , is: 

/ t ( p , Zi) = K ■ K + Ei ■ E“ + K ■ E? + K ■ E? , (7) 

where the spatial arguments of the fields have been suppressed and the * notation 
indicates that the complex conjugate is to be taken. The first term is constant at 
fixed Zd- The second term is necessarily non-linear in the susceptibility. Under weak 
scattering conditions, it is negligible compared to the third and fourth terms and will 
be disregarded hereafter. However, if the reference field is evanescent the scattering 
term becomes more important and eventually dominates as z d increases. The third 
arid fourth terms are conjugates that carry image information, which together will be 
referred to as the holographic intensity: 

1 (p, Zd) = E‘ y ■ E T a * + E T a ■ E£ . (8) 

It should be noted that I(p, Zd) may be positive or negative. 

Determination of r] is facilitated by taking the 2D (transverse) Fourier transform 
of the holographic intensity I(p,Zd) y i.e., 

55 I(P> z d) exp[zQ • p) . (9) 
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Substituting Eqs. (5) and (6) into (8) yields 


/(Q, z d ) = [//(Q,^>rXQ 4- q i - q r ) + £P(-Q, z d )fj *{- Q + q* - q r )] . (10) 

where rj* represents complex conjugate of the Fourier transform of rj. 

The expression for H{ Q, z d ) is: 

H(Q,z„) = g e? 4^-Z^ eX p{ife(q r -Q)-^(q,.)W}. (11) 

The tensor T 7 5 contains the integration over z' and is given in the appendix. It, is 
worthwhile to note the following properties of Eq. (10): 

1. When |q r — Q| > k n or | q r | > k 0 , the value of H (Q, z d ) decreases exponentially 
with z d . Similarly, when |q r + Q| > k 0 or |q. r | > k 0 , the value of H*(—Q,z d ) 
decreases exponentially with z d . 

2. The arguments of both rj and fj* are shifted by q, — q r . This allows the possibility 
for some of the super-oscillatory components of rj(Q) with |Q| > k 0 to be carried 
into the far field via homogenous waves. 

At this point one may generalize the result of Ref. 17, where the scattered fields 
in T1RM experiments are related to object structure, to include the effect of the 
holographic measurement. To understand the results it is helpful first to consider a 
heuristic argument. The objective is to determine the maximum value of |Q| such that 
the value of 77 ( Q) may be determined from measurements in the far-zonc. Consider 
only the first term in Eq. (10). In order to avoid exponentially decreasing H(Q ,z d ) 
with Zd . q r and Q must satisfy the conditions |Q — q, | < ko and |q r | < A: tJ . Recalling 
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that for the TIRM and PSTM imaging modalities |q ; | < nka (where n is the prism 
index of refraction), and setting Q - q r = k-ox where x is the unit vector in the x 
direction; it may be seen that |Q -j- q* — q r | is maximized when q, = nk^k. Thus, 
Q + — q, = (n + 1)&oX. Since the x direction was chosen arbitrarily, x may be 

replaced by any unit vector in the transverse plane. Therefore, in principle, f/(Q) 
may be determined when |Q| < (rtH- 1 )ko- This is enough information to construct 
a low-pass filtered version of the object. For comparison, the Fourier transform of 
the expression given by So 12 for the point spread function of a confocal microscope 
leads to a bandlimit of [Q[ < NA fc 0 , where NA is the numerical aperture (recall that 
NA < n) . 

3. Statistical model 

Real measurements are made on a space-limited array. For simplicity, assume that 
the detection surface has an area of L x L. lies in the plane 2 = z^. and that it has 
M x M pixels that give the total intensity at points {p mn }- Since the H functions 
are proportional to exp(— |Q|^d) when |Q| >> q r and |Q| >> q*, the intensity in the 
detector plane is effectively band-limited. Therefore Shannon’s theorem applies and 
the intensity in the detector plane can be represented by a discrete set of uniformly 
spaced samples. 18 The pixels are assumed to be small enough so that the intensity 
/ T (pi Zd) in Eq. (7) is essentially constant throughout each pixel (i.e., the detector 
does not introduce significant low-pass spatial filtering). It is also assumed that the 
detector has enough surface area so that nearly all of the scattered light is collected. 
Henceforth. M is taken to be an odd integer that is large enough so that aliasing is 
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negligible. This leads to the following discretization of the spatial frequency variable 
Q. The discrete spatial frequency is defined as Q mn = (/? m , /? n ), where 6 n is 

0 2tf o ^ M - 1 
On = , 0 < n < — - — 

2tt . , M - 1 

= — [n - M] , — — 4- 1 < n < M ~ 1 . (12) 

Note that 0*f_ n = — 3 n . 

In this discussion the semiclassical model for photoelectric detection is used, i.e., 
the electromagnetic field incident on the detector is treated classically, but the in- 
teraction with the detector is quantized in units of photoevents. 19 The intensities 
I T (p, z d ) and /(p, z d ) are considered to be deterministic quantities. The actual num- 
ber of photons counted in the mn th pixel (centered at transverse coordinate p Tnn ) 
is represented by the random variable N mn . The expected number of photons is 
N mn — WI T (fJ mn ,z d ), where W = (AaAt)hcfco, and the overbar is the expecta- 
tion operator, A a, A l, and ti are the effective pixel area (including the quantum 
efficiency), the measurement time, and Planck’s constant, divided by 27 t, respectively. 

With this scaling, Nq = WE r a - E r *, which is the expected number of reference 
beam counts from the first term in Eq. (7). It is convenient to define the random 
variable c mn = N mn - Nq, and thus cTif = WI mn , where l rnn is the holographic 
intensity from Eq. (8). 

It is clear from Eq. (10) that unambiguous determination of fj without the advan- 
tage of prior information (other than finite support) requires at least two measure- 
ments with different experimental parameters. Let those two measurements be de- 
noted (ci and C2). The estimate of i) is represented by 17, which is to be determined 
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by the simultaneous solution of the Eq. (10) for the two data sets with J(Q mn ) on the 
left hand side replaced by a measurement c. au jW . Due to noise, fj ^ f). Regardless of 
the probability density function under consideration, 7 mn = c mn /W , and Eq. (10) for 


the case of the two measurements takes the form: 


l Cl J 77 m 


V ftL ) 


/ 


= WHi(Q mn ) 


m mn T Q-il Rrl) 


\ 


\ ^?(Qmn d - 0.7 1 Qrl) J 


+ W H*{ Qm ) 


V m,M— » "b q»i Qr 1 ) 




(13) 
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(^2 )p+m,g+n 
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V fe) 


— W HliQp+m.q+n) 


2/p+m,(l+n J 


\ V* (Qm— m,Af— n + Chi ~ Rrl) J 

\ 

ij(Qp+m.()+n "b Qt2 — Qr2) 

^ P(Q p+m,q+n + Rt2 Qr'i) J 

t 


+ wh;( q p+m 

- m,q+M-n ) 


'7 (Qp+M — m,q+M—n + — q?- 2 ) 

^ V (Q p+M-m,q+M -n ~b Qj2 — Or 2 ) J 


where the integers (p + M — m) and (q -b M — n) are to be calculated modulo 
M. The integers p and q are chosen so that the arguments of f) and fj* in Eqs. 
(13) arc equal to within the discretization of $ given in Eq. (12). Thus, (p, q) « 


(L/27r) [qa — q r i — (q$ 2 — q r 2 )]- Note that in the case of PSTM q* — q r and p = q = 
0. In order simplify the notation, the equations (13) will be rewritten as follows, with 
the obvious symbol replacements. 


/ \ 

Cl 

= WH? 

f \ 
v+ 

+ WH{* 

( \ 
c 

V J 

( s 

C 2 

= W II 2 

K v + J 
( . > 

-5- WH.J* 

/ \ 
V*- 



p + , 


\v- 
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A. Gaussian PDF 


In this section the probability density function (PDF) of the inverse discrete Fourier 
transform (IDFT) of the spatial distribution of photocounts measured on the detector 
is determined. Consider a Gaussian noise model in which each value of N mn is statis- 
tically independent of the others and has a variance of o 1 . The equations are made 
more transparent by defining the zero-mean random variable b 7rLTl = N mn — Wl Tmn = 
c. mn — W I mn . It is convenient to represent all of the values of b mn as a single vector 
b = which will be called the zero-mean photocount vector. 

Since the photocounts in each pixel are statistically independent of each other, 
the PDF of b is: 

i r _i M-1A/-1 •) 

Ps(b) = (15) 

The vector b rnay be written in terms of its inverse DFT b: 

~ ( km Znl'| 

= S g “P {-#« [ M + m\ \ W 

Following the standard procedure for changing variables, 19 Eq. (15) may be used to 
compute the PDF as a function of b. Substituting Eq. (16) into Eq. (15) gives a 
sextuple sum. Making use of the identity: 

Y exp | —i2n T7 (k + p) } = MS(k - (M - p)) , (17) 

m=0 ^ M J 

where S represents the Kronecker 5. Eq. (15) becomes: 

— Af 2 M-1M-1 

l^B(b) = -srr E £ ^mn^M T Const . (18) 

17 n=0 m=0 

The PDF in Eq. (18) may be simplified. Since b is real and the n th and and [M — n) th 
points arc conjugate points in the M-point DFT, one it is seen that = Km- 
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This relation allows one to reduce the number of terms in the sum in Eq. (18) by 
almost half when it is written in terms of a minimal parameter set as follows: 


lnP-gCb) = 


-M 2 

2cr 2 


M-l 

2 


M - 1 


M — l 


M5 o + 2 £ {b Om Km + &m0 Ka) + 2 V 

m—l ■«=! m=l 


L* 


const . 

(19) 


Thus, each Fourier component included in this minimal parameter set is statistically 
independent of the others. If the variance of N mn were not constant with m and n, 
the various Fourier components would be coupled in the PDF. 

When there is more than one measurement of the photocount vector, presumably 
each measurement is statistically independent and therefore the joint PDF for all 
measurements is simply the products of the individual ones. If tq is the photocount 
vector of one measurement and b 2 is that of another measurement, the joint PDF for 
the two measurements is given by: 


PBB(b] , b 2 ) 

= ^b(bi)Ffe(b 2 ) and 

(20) 

^ 2 ) 

= P*$>. 

(21) 


Eqs. (14) and (19) may be used to characterize quantitatively how much infor- 
mation is may be obtained as a function of spatial frequency Q (in the object space), 
detection plane distance Zd, incident beam transverse wave vector q,, incident field 
e*, reference beam transverse wave vector q r and reference field vector e r . 

Since b mn — Cm n - WI mn , Eqs. (14) and Eq. (21) yield 

—W 2 M 2 ( 

In P mn (fj....,f}l,fi*_,r} + ) = — ||//i + (»} + “ r) + ) + *(r£ - ff*_)\ 2 

-r\H}{r)+ ~ fj+) + - f?-)| 2 } + const . (22) 
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In Eq. (22) r)_, ,f/ + are treated as independent variables. Eq. (22) is an example 

of the more general normal probability distribution InP(x) = const - (l/2)x F C _1 x 
where C is the covariance matrix of the random vector variable x and the superscript 
H is the Hermitian conjugate operator. The inverse covariance matrix in Eq. (22) is 
seen to be: 

0 0 £ s 

0 0 £* <r 

q e* 0 0 

£ c 0 0 

where g = iJftff* + e = H+H{ + and c = JJ+H+* + Eq. 

(23) has a symmetric block anti-diagonal structure. Its inverse is given by: 

0 0 c -e 

0 0 -e* g 

<L -£* 0 0 

— E Q 0 0 

The variance of the estimated quantities. |?/+ — ? 7 +| 2 and \f}_ — ?)_ | 2 , is given by the 
(2,4) and (1,3) elements, respectively, of the covariance matrix in Eq. (24). Since 
r]+{Q) = r^(— Q), only one of these is non-redundant. Thus the variance of the 
estimate is given by, 

^ = W 2 M-(e<, - ee*) ' (25) 

In addition, the covariances of the real and imaginary parts of fj + and ?)_ may be 
determined from Eq. (24) . Since the diagonal elements of the covariance matrix (24) 
are zero, the covariance of the real and imaginary parts are equal, and the real and 
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imaginary parts are uncor related. Thus the covariances of the real and imaginary 
parts of fj + and f/_ are equal to 1/2 of the values on the right-hand side of Eq. (25). 
It is interesting to note that according to Eq. (24) f)+ and i)_ are correlated. 


D. Poisson distribution case 

Although determining the PDF of the estimators 77 + , rf_ . fj + [as in equqatiori (22)] 
in an analytically tractable way for the Poisson statistics case is challenging, determin- 
ing the covariance matrix of the estimators as in Eq. (25) is relatively straightforward. 
Using the fact that (N mn — N mn ) 2 = N mn for a Poisson distributed random variable 
N mn , one may show that: 


C\C\ — CiCj CiC'2 — C 1 C 2 


c 2 cf - C 2 C 1 c 2 c* 2 - c 2 c 2 


U-4 r'V-1 \p.Vf-l a: 

2_> m =0 0 iV l,m 


M- 1 Eto 1 E" o 1 N- 


2 ,mn 


(26) 


The equality associated with the (1.1) and (2,2) elements in Eq. (26) is essentially 
identical to the result given by Goodman 13 in his discussion of stellar speckle inter- 
ferometry (the M -4 factor being due to the difference between the 1DFT and the 


DFT). 

One may use Eq.s (14) to express the variance \rj + - rj + \ 2 = fj + fj% — 77 + 77 + in terms 
of the quantities on the left hand side of Eq. (26). The result is: 


v+v+ - v+vl = 


1 


W 2 (q<; — 65*) 

/ \ T.T - |2 M - 1 M —\ 

n 21 EE 4 


M 4 


I it - 12 M — 1 M— 1 
1 1 1 


m— 0 n— 0 


M* 


(27) 


m=Q n=0 
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In order to bring Eqs. (25) and (27) into agreement, one must have No = a 2 , and 


Em E« A r i ,m» = Em En N 2>m n = M 2 N 0f which will always be true when I Tmn » I mn 
(i.e., the object creates a weak hologram) and will sometimes be true otherwise. Then, 
Eq. (27) becomes: 


fi+vX - v+vl = 


pN o 

W 2 M 2 {& - ee*) ‘ 


(28) 


C. Arbitrary noise case 


It is instructive to consider a measure of the expected error independent of the par- 


ticular noise model, so long as the noise is independently distributed (as was assumed 
in both cases above). One such measure is the squared I 2 norm of the system matrix 
inverse implied in Eq. (14). 20 That is, consider the matrix 


Ht Hr 
h} Hr 


(29) 


and the vector 


( . \ 

{ V- - V*- ) 


(30) 


related to the data through Eq. (14), f) = ,4b. The error in the estimate is necessarily 
bounded by the error in the data and the norm of the matrix inverse. Explicitly 


II A?>|| < U- 1 Ab 


(31) 


Consequently, the more singular A , i.e. the more poorly conditioned the system, the 
larger the magnification of the error in the estimate is expected to be. It is thus useful 
to consider ||vl -1 || 2 as an upper bound on the normalized error in the estimate. 
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4. Numerical results 


A standard mathematics package (MATLAB) was used to evaluate numerically Eq. 
(25) for PSTM and TIRM. In the examples shown here the variance of the esti- 
mate was calculated as a function of spatial frequency at five heights z d = (0.03 ; 
0.1, 0.3, 1.0, 10.0) A (where A — 27 t//cq). The measurements at greatest height be- 
ing effectively of the TIRM type rather than PSTM. In PSTM the reference field 
is the evanescent incident field, so that q 7 . — q* and e r — eh For the case that the 
field is measured at, a detector height, of 10A, it is assumed that the reference field is 
still of the same polarization and transverse wave vector as the incident field. Two 
independent polarizations of the incident- held are considered, TE (e l ;| k(q*) x z) 
and TM (e* || k(q») x k(q*) x z) modes with the following parameters: n — 2.0, 
q tl = (— 2.0/f-y, 0), and Az — 10 _4 A. For one set of examples q* 2 = (2.0fc O; 0), while 
for the second q ?2 — (0, 2.0&o). Thus the reference beams are anti-parallel for the first 
and perpendicular for the second. When q. (2 = (2.0A: 0 , 0) the experiment achieves the 
maximum possible transverse resolution possible (in the 2 -direction) . 

Figures 1 and 2 show the variance of the estimate fj( Q) as a function of spatial 
frequency. These “variance images” were normalized so that the smallest variance in 
any image is unity. The normalized variance images had a threshold set at 10 4 (i.e. 
values greater than 10 4 were set to that value) to reduce the dynamic range required 
of the display. Note that a value of 10 4 corresponds to an RMS error of 100 times 
that of the best determined value. 

It may be seen from Figure 3 that the squared norm of the matrix inverse is a good 
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indicator of which regions of the data space will be practically accessible. However it 
should be noted that this upper bound does give quantitatively different results from 
those obtained for the variance with particular noise models, and so || ^4 — 1 1| 2 cannot 
simply be taken as a substitute for the variance, nor can the noise model be ignored 
entirely. 

It can been seen in Figures 1. 2 and 3 that at Zd — A the variance and the norm 
of the matrix inverse increase quickly above the far-ficld cutoff |q* - Q| < fc 0 . The 
result indicates that the so-called far zone (the region where the evanescent fields 
are negligible) begins in the neighborhood of a single wavelength from the sample. 
If one hopes to take practical advantage of the high spatial frequency content of the 
near-field, measurements of the scattered field must be made within one wavelength 
of the sample. It should be noted that the shift of the Ewald sphere (or circle in this 
case) of reflection by the wave vector of the incident field remains evident even in the 
far zone. This shift may be used to advantage to circumvent the usual Abbe- Rayleigh 
limit when multiple incident fields of different wave vectors are used to sweep out a 
region of the Fourier space of the object outside the usual limiting circle of radius 
2A: 0 - 

The case of an incident wave that is less evanescent, q, — (— l.lfc 0 .0), was also 
considered, though figures are not shown here. The results were qualitatively similar 
to the above case and led to similar conclusions. That is, at Zd = A the variance of the 
estimate increases quickly above the far-fleld cutoff [in this case |Q| < (1.1 + 1 )k 0 ]. 
Of course, in such a case the incident wave does not provide much more resolution 
than a homogeneous wave with |Q| = fc 0 . 
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5. Conclusion 


An analysis of achievable resolution for optical near-field microscopy with has been 
presented. An error in the estimate of the computed dielectric susceptibility was 
obtained for Gaussian and Poisson noise models and a model independent measure 
of expected error was similarly obtained. The formulae presented here place limits on 
the spatial resolution that a given experiment may be expected to achieve for thin 
(or at least invariant along the z-axis), weakly scattering samples. 

The error in the estimate of the object structure was considered for measurements 
made at different heights above the sample. The region of the Fourier space of the 
sample practically accessible from these measurements is seen to decrease in size as 
the measurement plane is withdrawn to the far-zone. It was shown that the likely 
susceptibility estimate error increases dramatically for spatial frequencies above the 
far-ficld cutoff (|q, — Q| — &o), when the probe is even one wavelength away from 
the sample, indicating that the field does not contain much more information than 
is available in the far- field. The shift in the Ewald circle of reflection by the. wave 
vector of the incident field was seen to be manifest even in the far zone indicating 
that superresol ved imaging based on measurements taken in the far-zone alone is 
practically possible so long as the incident field is evanescent. 

Objects exhibiting complicated structure, not separable in the transverse and 
longitudinal coordinates, are also of great interest. In those cases, one may employ 
similar methodology and make use of the inverse scattering formalism described in 
Ref. 13 16 in combination with a noise model to determine the likely errors in the 
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object reconstruction. This is the subject of a forthcoming paper. 
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Appendix 


Expressions for the 1/2-space Green’s tensor may be found in the literature. 1 ’ 16 The 
tensor V 0 .q in Eq. (11) is obtained by integration of the Green’s tensor over the 
longitudinal coordinate of the object (z') and is given by the expression: 


r(q) = 



& XX +1l h VV ^xQyi^xx h-mi ) 

<h<ly( h XX hyy ) il h XX + Qx^yy 


\ 

|q| Qx^'xz 

|q|Qy/l IZ 


|q| <i*h zx |q|^ zx |q| 2 * zz J 


(Al) 


where q = (q x , q y ). The h functions are the result of multitplying the g functions in 
Ref. 16 by exp[i(k z (qi) — k; z (q))z'} and integrating from 0 to Az. Explicitly, 


hxx — ' ($i + -^2^2) 

Kq 

= dsMs Ifr-w 

Kq 

llyy = (flj+JeA) 

hzi = + R2f)2) 

Kq 

In E 

h zz = ^(0: - R 2 O 2 ) , (A2) 
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where Ri(q) and R 2 (q) are the reflection coefficients given by 


R i(q) = 


fc,(q) - k' z ( q) 
k z ( q) + ^(q)’ 


(A3) 


and 


„ , , _ k 'M) - »fc«(q) 

2(q) k’(q) + nk z (qy 


with k' z { q) = yjn 2 k 2 — \ q| 2 . The functions 9 ^ and 9 2 are given by: 


0 i(q) = i 


.1 - exp (i[k z (qi) - /c*(q)]Az) 


Mq*) - Mq) 


, . .1 - exp(i[k z {qi) + fc z (q)]A 2 ) 

02(q) = * tM+k^) 


(A4) 


(A5) 
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Fig. 1. The base- 10 logarithm of normalized PSTM susceptibility Fourier space 
variance for detection heights of (0.03, 0.1, 0.3, 1.0, 10.0) A as labeled above. The 
coordinates are spatial frequency in units of ko = 27r/A. The data are taken 
to consist of two scans with incident wave vectors of ±2/c 0 x, where x is a unit 
vector. Results of illumination by TM modes are shown on the left and TE on 
the right. The plots are normalized so that the minimum variance of any point 
is 1.0. The linear color scale runs from 0 to 4 and is shown at the bottom of the 
figure. Note that the range on each figure is different. Values of the normalized 
variance greater than 10 4 were set equal to 10 4 for clarity in display. It can be 
seen that at z d = 1.0A the effective information content is already similar to 
that of the far-ficld limit, in which only the homogeneous modes are detected. 

Fig. 2. Same as Figure 1, except that the data are assumed to consist of two 
scans with incident wave vectors of —2 fc c x and 2/c 0 y, where y is a unit vector. 

Fig. 3. The logarithm of the squared I 2 norm of the system matrix inverse for 
the cases of measurements made in the two closer planes. The four top plots 
show results for the case of counter propagating incident evanescent waves, 
while the bottom four plots show results for the case of the same orthogonal 
incident wave vectors discussed above. Coordinates are again spatial frequency 
in units of ko- The linear color scale running from 0 to 4 is shown at the bottom 
of the figure. 
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